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Abstract 



. The potential energy curves, permanent and transition dipole moments, and the static dipo- 



I lar polarizabihty, of molecular ions composed of one alkali-metal atom and a Strontium ion are 

Ph. 

determined with a quantum chemistry approach. The molecular ions are treated as effective two- 

^ • electron systems and are treated using effective core potentials including core polarization, large 

^3 ■ 

gaussian basis sets, and full configuration interaction. In the perspective of upcoming experiments 
aiming at merging cold atom and cold ion traps, possible paths for radiative charge exchange, 

> 

■ photoassociation of a cold Lithium or Rubidium atom and a Strontium ion are discussed, as well 

<N , 

. as the formation of stable molecular ions. 
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I. INTRODUCTION 

As illustrated by two recent review papers , the success story of researches on ul- 
tracold neutral molecules is developing for more than ten years since their first detection in 
an ultracold cesium gas js-Sj. Ultracold molecules indeed offer exciting prospects for funda- 
mental quantum physics and physical chemistry, as well as for interdisciplinary physics and 
possible applications and for technology. The ultimate control of their external degrees of 
freedom down to kinetic energies equivalent to temperatures far below 1 mK, and of their 
internal degrees of freedom at the single (lowest) quantum level is a crucial path 

towards the realization of a degenerate quantum gas of ultracold molecules in their absolute 
ground state. Ensembles of ultracold molecules in well-defined quantum states allow for the 



observation of cold collisions between atoms and molecules 
tary chemical reactions is now at reach 
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. The control of elemen- 



19l-l25|. while new light on the quantum behavior of 



few-body systems is brought through the demonstration of predicted universal laws |26l. |27]. 
Many-body phy=ic= could be engineered wjth ult.acold pola. molecule, serving a. n.odel. 

for spin lattices with tunable interactions [28], or as quantum bits in the perspective of the 



long-term achievement of a quantum information device [29|, |30| . The unprecedented accu- 
racy on transition frequency measurements induced by their ultra-low kinetic energy makes 
polar molecules suitable for testing fundamental theories, like for instance the variance of 
the electron-to-proton mass ratio rrie/mp 3lN33| or electron dipole moment (EDM) 171. |34|. 
Besides, molecular ions are well-known for their crucial role in the evolution of many 
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36| , and it is highly desirable 



environments like stellar atmospheres or interstellar media 
to create samples of cold molecular ions in order to learn all the details of their quantum 
dynamics and to control their reactivity 
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. In this matter, electromagnetic traps for 
ions have long been used |37|. Novel developments occurred with the demonstration of 
sympathetic cooling of molecular ions with laser-cooled atomic ions thus formin g so -called 



ion crystals where individual ions are kept at a well-defined positions in space l38l|. It is 
nowadays possible to drive such ions into selected rotational level through additional laser- 
pumping mechanisms 
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4l| . Just like neutral molecules, cold molecular ions are promising 



candidates for rrie/mp variance j42j or EDM 43|, measurement. 



For such studies, the prerequisite step is obviously the formation of a cold and trapped 
sample of the desired species of molecular ions which would be sympathetically cooled by 
laser-cooled and trapped atomic ions. It is now established that excited ions in the trap 
and surrounding molecular species introduced on purpose easily react together to produce 



molecular ions, as demonstrated for the first time in Ref. 



45l | with Mg"*" and H2 yielding cold 



and trapped MgH"*" ions. With the developing experimental trend of merging traps of laser 



cooled atoms and ions like Yb-Yb"*" 



46j, Rb-Yb- 



47|, 



48^, Rb-Ba+ 



49], new opportunities 



could be investigated to create diatomic molecular ions in situ via e.g. the photoassociation- 
stabilization sequence A + + hu ^ (AB^)* — > AB^ + hu', or via magnetoassociation 

n n n n 

involving Feshbach resonances [50|, just like routinely done with cold atoms [2|, |51l, l52 |. 
Charge exchange between the neutral and ionic atomic species should also be considered, as 
explored in the Na-Ca"*" system 0, 

The main drawback for accurate theoretical modeling of these formation mechanisms 
often arises from the poor knowledge of the molecular structure of such mixed species beyond 



their electronic ground state. As mentioned in our previous work [5J] (hereafter referred to as 
paper I), this is particularly true for the diatomic systems, either neutral or charged, formed 
from an alkali-metal atom and Strontium. The Strontium ion is indeed an attractive system 
for trapping and laser-cooling experiments (see for instance Refs. 55|, |56|). Following paper 



I, the present work is devoted to the theoretical determination of the electronic structure of 
(y4-'^''Sr)+ ions (where A= ^Li, ^^Na, ^^K, ^^Rb, ^^^Cs) using the method briefly recalled in 
Section ini Results for potential energy curves and permanent and transition dipole moments 
are given in Section UTT] and HVl respectively. In Section |V] we explore the possibility to create 
cold (A-Sr)"*" ions via photoassociation followed by spontaneous decay, and the expected 
reactivity of these ions with surrounding cold neutral atoms. Atomic units for distances 
(1 a.u.= aO = 0.052 917 720 859 nm), energies (1 a.u.= 2R^ = 219474.63137054 cm"^), 
and dipole moment (1 a.u. =2.541 58 Debye) will be used throughout the paper, except 
otherwise stated. 



II. COMPUTATIONAL APPROACH 

The present calculations actually represent an intermediate step of those performed in 
paper I for neutral A-Si molecules, which were modeled as effective three-electron systems. 
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The A-Sr"*" ion has two valence electrons, and its electronic structure is, at least in terms of 
values of quantum numbers, similar to the one of alkali dimers. Therefore we describe the 
A-Sr"*" ion as an effective two-electron system following the same steps than in paper I and 



in Ref . 



571], which are briefly recalled below. 



Each of the closed-shell ionic cores A'^ and Sr^"*" are represented by an effective core po- 



tential (ECP) from Durand and Barthelat |58|, l59j and Fuentealba et al. [60|, l6l| respectively. 



completed by the core-polarization potential (CPP) depending on the angular momentum 



62| and revisited by Guerout et al. 



63|. 



^ of the valence electron, along the lines of ref. 
The valence electron wave functions are expanded on a Gaussian orbitals centered at each 
ionic core. The molecular orbitals are determined by restricted Hartree-Fock single electron 
calculations, including the CPP term yielding the potential curves for the relevant molec- 
ular double-charged cations. Then a full configuration interaction (FCI) is performed for 
each relevant molecular symmetry, providing potential curves, and permanent and transi- 
tion dipole moments. The calculations are performed through an automatic procedure js^l 
based on the the CIPSI package (Configurationjnteraction by Perturbation of a multicon- 
figuration wave function Selected Iteratively) [64!]. As mentioned in paper I we use for the 
Sr+ ion a large uncontracted basis set {7sQp7d} derived from the {5s5p6(il/} used in Ref. 



65|, and {8s7p4(i}, {7s6p5(i}, {7s5p7(i}, {7s4p5(i}, and {7sApbd} uncontracted basis sets 



for Li, Na, K, Rb, and Cs, respectively. All exponents of the Gaussian orbitals, as well as 



the cut-off radii and the core polarizabilities |66|, |67| appearing in the CPP term for the 
alkali species, were already reported in paper I and are not repeated here. The cut-off radii 
for the Sr^"*' CPP related to the {IsQpld} basis set are slightly adjusted compared to paper 
I, i.e. ps = 2.13035, ps = 2.183, ps = 1.70616, in atomic units. The quality of the basis sets 
is settled through several criteria. First of all, the cut-off radii are adjusted to reproduce the 
experimental energy of the lowest levels of the alkali atoms and and of the Sr"*" ion (see Table 
3 in paper I). Next, as the lowest dissociation limit of the A-Sr"*" ions is A++Sr, the lowest 
energy level of the -effective two-electron- Sr atom calculated by the FCI above should be 
well reproduced. As shown in paper I, the present basis set is quite efficient in this respect, 
as the Sr ground state and 5s5p^P, and 5s5p^P excited state energies are calculated with 
a difference of 185 cm~^, 269 cm~^, and 434 cm~^ compared to the experimental values, 
respectively. The corresponding excitation energies of the 5s5p^P and 5s5p^P levels are 
obtained with a difference of 87 cm~^, and 262 cm~^ compared to the experimental values. 



which represents an accuracy comparable to the one achieved in the recent high-accuracy 



calculations on Calcium by Bouissou et al. 



68|. 



In addition, the quality of the electronic wave functions of the Sr"*" ion can be checked 
against the calculated static dipole polarizability, which has been studied theoretically for 
several energy levels in various contexts: the broadening of spectral lines due to collisions 
with hydrogen atoms 69|] , the long-range interaction between the Sr+ ion and various atomic 



species 



70, 



7l|, or the search for stable and reproducible atomic frequency standards based 



on the 5s — 4^5/2 clock transition in ^^Sr"*" 72l-l74l| . All these calculations are based on the 
summation over oscillator strengths at various levels of approximations, namely including 

ItiI or not 



transitions towards the continuum [70|, 171| or not |69|, |73|, 17^. Here we use the finite field 
method 75| for the contribution of the valence electron, as in Ref. 76| and in our previous 



work on alkali atoms [77[. More precisely we calculated the energies of the 21 + 1 components 
of an atomic state with angular momentum ^ in the presence of an electric field varying from 
0.002 to 0.004 a.u., and extracted the polarizability from the quadratic dependence of the 
averaged values against the electric field magnitude. Finally we add the contribution of the 
Sr^+ core (5.67 a.u after Ref. [ggI), as pointed out in Ref. 78|. Our results are reported in 



Table [H for the 5s, 4(i, 6s, 5(i, and 5p levels. Our va 
with all other values except the oldest one of Ref. 



the most recent ones of Ref. 



ues for the 5s and Ad states agree well 



76| , and are in reasonable agreement with 



7l| for the excited p and s levels. The underestimation visible 



in our calculations for the d levels could be due to the absence of / orbitals in our basis sets. 

We calculated with the same approach the polarizability of Sr in its ground state (Table 
H]) which agrees with less than a 10% difference with the experimental value of Ref. |79| . as 



well as with the recent theoretical value of Ref. [80| . All other theoretical determinations are 
larger by 20% or more. A new experimental determination of electric dipole polarizability 
of strontium then seems desirable to discriminate among all determinations. 



III. POTENTIAL CURVE CALCULATIONS 

An overview of the results obtained from the procedure above for the potential curves 
of the A-Sr"*" ions is provided in Figs. |2] to El while the main spectroscopic constants are 
reported in Tables HTl to IVTl In the tables, we recall for convenience the dissociation energies 
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State 


This work 


[71] 


[73J 


m 


m 


[69] 


[76j 


5s 


92.7 


89.8 


84.6 


91.3 


91.47 


93.1 


130.5 


Ad 


53.3 


61.77 


53.7 


62 




57 




5p 


-49.2 


-23.13 








-32 




6s 


1305 


1089 












5d 


1294 


2099 












6p 


-2016 


-2056 














This work 


[79] 


[81J-DF-HF 


[82j 


[83] 


[76] 


[80] 




200 


186±14 


228-241 


186 


225 


245.7 


197.2 



TABLE I. Static dipole polarizabihty of 5s Ad, 5d and 6p states of Sr'^ ion (in atomic units) 



compared with other calculated values. We added the Sr-^"*" polarizabihty of 5.67 a 



.u. leel] 



to the 



values of Refs. The indexes DF and HF hold for Dirac-Fock and Hartree-Fock level of 

calculation, respectively. The static dipole polarizabihty for the Strontium ground state is also 
displayed, compared to several theoretical values, and an experimental result [79|. 



of the asymptotic limits as displayed in paper I. Let us note that no other results for these 
systems have been previously published. The corresponding numerical data are provided as 
additional material accompanying the paper. 

For all ions, the ground state X^S+ of the system correlates at large distances to a 
ground-state neutral Sr atom and a closed-shell alkali-metal ion. The ordering of the well 
depths is the same than with neutral A-Si molecules the deepest wells are by far the 
LiSr+ and NaSr+ ones (11157 cm~^, and 8330 cm~^, respectively), while the well depths 
of KSr+, RbSr+, CsSr+ are significantly smaller (Fig. H]). This can be understood as the 
energy spacing between the two lowest dissociation limits increases with increasing alkali- 
metal atom mass: the X^S+ state is less and less repelled by the next (2)^S+ state, and 
the X state behaves more and more like a van der Waals molecule composed of two closed- 
shell systems A'^ and Sr(5s2). The equilibrium distances are comparable to the typical 



values of alkali dimers (see for instance Ref. 



841 ] where a sketch of potential curves for all 



heteronuclear alkali-metal dimers is displayed) as the molecules are formed by two open- 
shell partners. At the second dissociation limit, one valence electron is attached to each 
atomic partner, giving rise to a singlet and a triplet state at short internuclear distances. 




FIG. 1. Potential energy curves of the X S+ ground state of ^-Sr^ ions, all referring to the energy 
origin fixed at the A+Sr"*" dissociation. 



Despite this similarity with the lowest electronic states of ground state alkali-metal pairs, 
the situation is here strikingly different: the lowest triplet state (also referred to as 

the a^S+ state) potential curve is strongly attractive while the accompanying singlet state 
(2)^E+ (also referred to as the A^S+ state) potential curve is mostly repulsive, due to the 
strong repulsion induced by the X^S+ ground state. The (2)^S+ would actually represent 
a strong dissociation channel induced by spin-orbit coupling for the (l)^n state if collisions 
between alkali-metal ions and neutral Strontium were studied. However, as a reciprocal 
effect of the larger and larger energy spacing between dissociation limits invoked above, the 
heavier A-Sr"*" molecule, the deeper long-range potential well in the (2)^2"*" potential curve. 

The (l)^S'*' is significantly deeper (several thousands of cm~^) than the one of alkali- 



metal dimers (generally having a few hundreds of cm ^ depth, see Ref . [84| again). Indeed, 



the charge-induced dipole long-range attraction (varying as i?"^, where R is the internu- 
clear distance) between Sr"*" and the alkah-metal atom is stronger than in the alkah-metal 
pairs (varying as R~^), and dominates the (repulsive) exchange interaction down to short 
distances. As the static dipole polarizability (determining the strength of the long-range A- 
Sr"*" interaction) increases with increasing alkali-metal atom mass, the depth of the (l)^!!"*" 
increases, and becomes even larger than the X^S"*" well depth for KSr"*", RbSr"*", CsSr+. Note 



85|. 



that this trend is similar to what is also observed for the RbBa"*" ion 

The combination of the excitation energies of Sr+, of Sr, and of the alkali- metal atoms 
induces the main differences for excited molecular states among the five different A-Sr"*" 
species. For instance, the LiSr+ molecule is the only one where the dissociation limit with 
the alkali-metal atom excited in its lowest p level has an energy larger than the one involving 
Sr(5s5p^P). The Li(2p)+Sr+(5s) dissociation limit actually lies only 178 cm~^ above the 
Li(2s)+Sr^(4(i) one. Similarly, the LiSr"*" and NaSr^ species are the only ones in which 
the dissociation limit involving the Sr"'"(4(i) level has a binding energy smaller than the one 
involving Sr(5s5p^P). In this respect, the KSr"*", RbSr+, CsSr+ exhibit similar patterns in 
their excited states. Let us note that due to these differences, the lettered labels A, B, C, D, 
E, a, b, c of the molecular states, which is historically related to their energy ordering, will 
not always be attached to the same dissociation limits in the A-Sr"*" series. The consequences 
of such characteristics for the dynamics of the A-Sr"*" systems will be further addressed in 
Section |V] below. 

Among all the excited states, there are two cases where the corresponding potential curves 
exhibit remarkable features. First the (2)'^S+ curve is strongly repelled by the attractive 
(l)'^S''" one so that it strongly interact through an avoided crossing with the (3)^S+ curve 
located in the 9-12 a.u. distance range. For all species but LiSr+, it induces a strong mixing 
of the and states of Strontium, starting from large distances. This will undoubtedly 
lead to significant predissociation of the (3)'^S+ bound levels, without charge exchange. This 
creates a second long-range well in the (2)'^S+ curve, which almost disappears in CsSr+ due 
to the continuous decrease of the ionization potential of the alkali species with increasing 
mass. In the LiSr"*", this would actually represent a channel for charge exchange, as the 
(3)3S+ is correlated to the Li(2s)+Sr+(4(i) limit. 

Second, the pair of states (3)^S+ and (4)^S+ show a similar pattern while slightly more 
complicated. For the KSr"*", RbSr+, and CsSr+ species, these two states dissociate toward 
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the neighboring A'^+St{^D) and y4++Sr(^P) hmits, resuhing into an avoided crossing in the 
12-15 a.u. distance range, with a second long-range potential well in the (3)^X1"'" curve. In 
LiSr'^ a similar pattern is predicted between these two states, except that their interaction 
corresponds to an exchange of excitation for the Li atom, due to the different ordering of 
the excitation energies in the various species. In between, the NaSr"*" ion should exhibit 
a curious feature: the interaction between these two states results from the proximity of 
the Na(3s)+Sr^(4(i) and Na"^+Sr(-'^Z}) limits, so that two successive avoided crossings are 
visible. Our results suggest that these crossings are quite small, so that it could well be 
that these two states behave independently in a collision process. The (4)^S+ is expected to 
have a double- well structure as well, which disappears in RbSr"'' and CsSr+. These features 
are actually at the limit of the capabilities of our computations and would require a much 
larger basis set for an appropriate representation. 



IV. PERMANENT AND TRANSITION DIPOLE MOMENTS, STATIC DIPOLE 
POLARIZABILITIES 

Permanent dipole moments (PDM) are relevant to the present study, i.e. if molecular 
ions in the X^S+ ground state or in the a^S"*" metastable state are created in experiment 
manipulating cold Sr"*" and alkali-metal atoms: PDMs will determine the ability to perform 
vibrational or rotational cooling following for instance methods recently demonstrated for 
MgH"*" 40(1 and HD"*" 39| . The PDM is here calculated in the body-fixed frame with respect 
to the origin placed at the center-of-mass of the molecule, and is therefore characterized by a 
linear divergence at large internuclear distances expressing the increasing distance between 
the negative and positive center-of-charge. Our results are reported in FigJ7]for these two 
electronic states and for all A-Sr"^ species, confirming this pattern. The molecular axis 
is oriented along the Sr+-A direction, so that the positive sign of the ground state PDM 
indicates an excess of electron charge on the alkali-metal atom, which is less pronounced as 




R (Bohr radii) 

FIG. 2. Potential energy curves of the LiSr"*" molecule computed in the present work. full 
red lines; ^S^: dashed red lines; ^H: full blue lines; ^H: dashed blue lines. Vertical arrows suggest 
possible photoassociation (upwards) of cold Li and Sr+ particles and spontaneous decay paths 
(downwards) to create cold LiSr+ molecules. Note that curves of A symmetry are not drawn for 
clarity sake. 

its mass increases. This is compatible with the decrease of the electron affinity of alkali- 
metal species with their increasing mass. The PDM is almost linear in the region of the 
minimum of the ground state, so that radiative transitions among vibrational levels levels 
should be possible. Its magnitude is also significant, so that rotational transitions should be 
observed. In contrast the sign of the triplet PDM, while mostly positive, cannot be easily 
interpreted over the whole distance range. Note that its dissociation limit is different from 
the ground state one, so that its divergence at large distances is different. Its magnitude is 
almost constant (and even close to zero for CsSr+), which therefore is probably not a good 
candidate for cooling of internal degrees of freedom. As a complement to the description 
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FIG. 3. Potential energy curves of the NaSr"^ molecule computed in the present work. Line types 
as in Figj2j Note that curves of A symmetry are not drawn for clarity sake. 

of the properties of the A-Sr"*" molecules, we computed their static dipole polarizabilities 
(SDP), which will be relevant for instance when considering the long-range interaction with 
surrounding cold neutral alkali-metal atoms in a merged cold ion/atom trap experiment. 



As pointed out in Ref. 



78|, the SDP includes a core contribution (5.67 a.u. here, see Table 



II of paper I) in addition to the valence contribution, the latter being evaluated using the 
finite field method [tS]. Both parallel a\\ (along the internuclear axis), and perpendicular 
a±_ components (see FigJH]) ex^ 



dimers 



771 ] ■ alkali hydrides 



libit an i?-variation similar to the one obtained for alkali 

In all 



for MgH"*" 871], or for the hydrogen molecule [88 
these systems, a\\ has a maximum at a distance larger than the equilibrium distance Re of 
the system, while a± always has a smaller magnitude than ay, and monotonically increases 
towards the asymptotic limit. Around the equilibrium distance, a\\ ranges between 200 and 
250 a.u., and between 120 and 150 a.u.. As for the well depth the maximum of a\\ increases 
as the mass of the A-Sr"*" molecule decreases. At large distances, both components converge 
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FIG. 4. Same as Fig. [3] for KSr+. 



toward a unique limit corresponding to the sum of SDP of the Sr ground state and the (small) 
SDP of the alkali-metal ions. Finally, the electronic transition dipole moments (TDM) are 
relevant when considering for instance the possibility for radiative charge exchange (RCE) 



in these systems 



50 



or for photoassociation (PA) and spontaneous emission (SE) (see 
next section). We have computed all TDMs of all A-Sr+ species, whose numerical data are 
provided through the EPAPS system. In order to illustrate them, we chose the LiSr"*" and the 
RbSr+ molecules as being representative of the A-Sr"*" species (see Section [In|) . We represent 
in Figs. [9] and fTOl the TDMs involving X^'E'^, (2)^S+ and a'^S+ as the lowest states. For 
the sake of clarity, we will use both numbered (in parenthesis) and lettered (when available, 
in brackets) notations for the molecular states. Various patterns can be analyzed in these 
figures. 

The most obvious one is when the TDM function connects two molecular states which 
dissociate towards atomic states connected by a resonant transition, as for the (1)^S+[X] — 
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R (Bohr radii) 



FIG. 5. Same as Figi2]for RbSr+. 

(2) in[L)] and (l)iS+[X]-(4)iS+[£;] in RbSr+: the TDM function converges towards the Sr+ 
{^S P) atomic TDM. The i?- variation of the former suggests that the molecular states [X] 
and [D] keep a similar electronic character, while the abrupt variation of the latter reflects 
the strong configuration interaction of the [E] state, induced by the avoided crossing with 
the (3)^S^[C] state. It is actually clear on Fig JTOT a) that the [C] and [E] state exchange 
their electronic character at the position of their avoided crossing around 14 a.u.. The other 
displayed TDM functions for - (2)iS+[A], {iy^+[X]-{iyU[B], and 

(3) ^E+[C] vanishes at large distances, with strong variations at short internuclear distances 
again related to strong configuration interaction. The same analysis holds for transitions 
involving the [A] state. In particular, the states connected to [A], i.e. [C], [D], and [E] have 
vanishing TDMs at large distances, as they correspond asymptotically to different charge 
states of the separated atoms. A similar reasoning can be made for LiSr+, where the results 
reflect the different ordering of the asymptotic limits compared to the other species. For 
both species represented in [Hb and [TUb. the TDM functions for the (l)'^S+[a] — (2)'^S+[c] 
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FIG. 6. Same as Fig. [3] for CsSr"*". 



and (l)^S+[a] — (3)^S+[e] exchange their character at the distance of the avoided crossing 
between the [c] and [e] states around 9 a.u. and 11 a.u. for LiSr+ and RbSr+, respectively. 

V. DISCUSSION: RADIATIVE CHARGE EXCHANGE, PHOTOASSOCIATION, 
AND STABILITY OF A-SR+ MOLECULES 

The experiments aiming at merging a cold Sr+ trap and a cold alkali-metal atom trap 
involve molecular states connected to the second dissociation limit, i.e. the (2)^S+[A] and 
the (l)'^S+[a] states. Despite the strong variation of the [X] — [A\ TDM between 10 and 
20a. u., it is most probable that the RCE process during a cold collision between ground 
state alkali-metal atom A and Sr+ will be weak as the TDM rapidly vanishes with R. This 
is actually the trend already reported in Ref. jsS] for Na-Ca+ collisions. However, as it is 
possible to trap clouds of laser-cooled Sr+ ions for hours, it may be possible to observe the 
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FIG. 7. Permanent dipole moments (a.u.) of (a) the X^S+ ground state and (b) the a'^S"'" 
metastable state of the ^-Sr^" molecules computed in the present work. The origin is taken at 
the center-of-mass of the molecule, considering the ®''Sr, ^Li, and ^^Rb isotopes. The position of 
equilibrium distances are shown with vertical dashed lines. 

slow decreasing of Sr"*" ions, while the A'^ alkali-metal ions will manifest as dark spots in 
the ion crystal. The trapped ions are continuously illuminated by the cooling lasers, so 
that excited alkali-metal atom A* in their lowest p state are present. In this perspective, 
the Li-Sr"*" exhibits a peculiar situation: as noticed in Section [XTTl the Li(2p)+Sr+(5s) limit 
lies only 178 cm~^ above the Li(2s)+Sr+(4(i) one. The transition between the [D] and [E] 
^S"*" enhanced by their mutual avoided crossing may lead to a strong excitation exchange, 
leaving Sr"*" ions in their 4d state, insensitive to the cooling laser frequency. More generally, 
the asymptotes involving the resonantly excited Sr+(5p) lie outside of the energy range of 
the presented figures, where a large number of potential curves probably offer many channels 
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FIG. 8. (a) Parallel and (b) perpendicular components of the static dipole polarizability (in a.u., 
1 a.u. =0.148 184 709 3 A'^) of the X^S^" ground state of the A-Sr"*" molecules computed in the 
present work. The position of equilibrium distances are shown with vertical dashed lines. 

for charge exchange. It will be interesting to check the variation with the intensity of the 
cooling lasers of the ion number in such experiments, which would require an extension of 
the present treatment. 

Prospects for efficient PA of A-Sr+ pairs, and stabilization into cold and trapped ASt^ 
molecular ions look promising, according to Figs H] and El The Li-Sr"*" and Rb-Sr+ molecules 
are taken as representative examples, due to the above mentioned peculiarity of Li-Sr"*", 
and the use of Rb atoms in several merged ion trap/cold atom trap experiments. For 
instance, the excitation of the (2)^11 state from the entrance channel Li(2s)+Sr"'"(5s) with 
a wavelength slightly larger than the one of the 2s — 2P transition in Li could take place 
either from large distances, or from the inner turning point of the (2)^S+ potential curve. It 
could lead to an efficient decay down to the low-lying levels of the ground state. A similar 

16 




FIG. 9. Selected transition dipole moments between (a) singlet states, and (b) triplet states, for 
LiSr+. Molecular states are labeled with letters as in Table [III 



mechanism is predicted for the formation of metastable (1)^S+ ions in deeply-bound levels 
via the excitation of the (3)^S^. PA of Rb and Sr"*" pairs are expected using a wavelength 
detuned to the red of the — 5s5p^P Sr"*" resonance. An efficient decay is expected to take 
place from the inner turning point of the (4)^S+. Another path could be the excitation at 
the turning point of the (2)^S+ toward the (2)^11, which would be stabilized into the lowest 
ground state levels. No triplet states will ba accessible in such a scheme. But the spin-orbit 
coupling have been neglected so far. It is quite large in the Sr"*" ion (280.34 cm~^ for the 



4d level, and 801.46 cm ^ for the 5p level 89|), and comparable in some cases with the 



one of alkali-metal atoms (0.29, 17.2, 59.7, 237.9, and 551.3 cm~^, for Li, Na, K, Rb, and 
Cs, respectively). This will induce additional couplings between molecular states of various 
symmetries. For instance, the spin-orbit interaction will strongly mix the {2yU and the 
(3)^E+, probably leading to strong predissociation channels with charge exchange into the 
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FIG. 10. Selected transition dipole moments between (a) singlet states, and (b) triplet states, for 
RbSr^. Molecular states are labeled with letters as in Table IVl 

Rb++Sr(^P,'^ outgoing channels. A detailed description would require further quantum 
chemistry calculations which are beyond the scope of the present paper. 

Once y4Sr+ ions are formed, the important issue of their stability against charge ex- 
change with the surrounding cold alkali-metal atoms is raised, just like in the proposal of 



Ref. 90|] looking at possible sympathetic cooling of molecular ions by laser-cooled atoms. 
In our case, the ionization potential (IP) of alkali-metal atoms should be larger than the 
IP of the related neutral molecule. Combining the present results with those of paper I 
for the ground state of neutral ASi molecules, we estimate their IP by the energy dif- 
ference between the minimum De of their ground state potential curves according to: 
P(AST)=IP{A)+Df;{AST)-De{AST~^). Using the experimental IP for alkali-metal atoms 



89|, namely 43487, 41449, 35009, 33691, and 31046 cm^^ for Li, Na, K, Kb, and Cs re- 



spectively, we find IP (LiSr) =34848 cm'^, IP(NaSr)=34716 cm^^, IP(KSr)=31198 cm-\ 

18 



IP(RbSr)=31126 cm ^, IP(CsSr)=29416 cm ^. Therefore we predict that the ASr~^ ions 
will be all stable against charge exchange with alkali-metal atom A. 
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State 


Label 


Diss, limit 


-E'(i2ss(a-u.) 


(5(cm-i)) 


Re 


De 






X 


Li++Sr(5s2i5) 


-0.61548 


-184 


6.37 


11126 


246 


(2)iS+ 


A 


Li(2s)+Sr+(5s) 


-0.60350 





15.33 


277 


35.4 




a 






77 


7.48 


4865 


163 


(i)3n 


b 


Li++Sr(5s5p3p) 


-0.54890 


-274 


5.90 


11087 


268 


(2)3s+(*) 


c 


„ 






6.90 


-5086 


190 


2""^ well 




„ 




„ 


12.1 


-3205 


76 


barrier 




„ 


•)•} 


„ 


9.80 


-2652 




(i)in 


B 


Li(2s)+Sr+(4(i) 


-0.53641 




7.14 


6150 


175 


(3)^S+(*) 


C 




11 




7.32 


-5479 


130 


2"'^ well 




„ 


77 


„ 


13.8 


-1297 




barrier 




„ 


77 


„ 


12.9 


-1259 




(3)3S+ 




„ 


77 


„ 


9.54 


5722 


69 


\ J 




„ 


77 


„ 


7.21 


6098 


58 


(l)iA 

V / 




„ 


77 


„ 


6.92 


7851 


194 


(1)3A 






77 




6.97 


8003 


193 


(2)in 


D 


Li(2p)+Sr+(5s) 


-0.53559 




6.64 


2460 


198 


(4)iS+(*) 


E 




77 




7.87 


-1183 


260 


2"^^ well 






77 




12.59 


-1255 


157 


barrier 






77 




9.32 


-h591 





TABLE II. Main spectroscopic constants of the lowest electronic states (with their conventional 
spectroscopic label) of the LiSr+ molecule: the equilibrium distance Re (in atomic units), the 
well depth D^. (in cm~^) and the harmonic constant uje (in cm~^). For convenience, we recall the 
calculated dissociation energies EfUss, and their difference with the experimental ones [sj]. The (*) 
sign denotes double-well potentials, for which the barrier height and the energy at the minimum 
of each well are displayed (in italic) with respect to the dissociation limit, besides their position in 
R. 
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State 


Label 


Diss, limit 


Edissia^.n.) 


5(cm-i)) 


Re 








X 


Na++Sr(5s2) 


-0.61548 


-184 


6.92 


8330 


121 


(2)iS+ 


A 


Na(3s)+Sr+(5s) 


-0.594215 





14.8 


319 


21 




a 




77 


77 


8.03 


4363 


83 


(i)3n 


b 


Na++Sr(5s5ij3p) 


-0.54890 


-274 


6.58 


7491 


127 


(2)3s+(*) 


c 




77 


77 


7.43 


-2765 


98 


2"^^ well 






?7 


77 


13.1 


-2417 


39 


barrier 




77 


77 


77 


9.8 


-1017 




(2)3n 




Na++Sr(5s4d3D) 


-0.52991 


343 


7.80 


5542 


85 


(1)3A 




77 


77 


77 


7.45 


7091 


98 


(3)3S+ 




77 


77 


77 


9.90 


4297 


148 


(i)in 


B 


Na(3s)+Sr+(4d) 


-0.52712 





7.64 


5911 


90 


(3)iS+(*) 


C 


77 


77 


77 


7.95 


-5402 


86 


2"^^ well 




77 


77 


77 


13.4 


-2536 




barrier 




77 


77 


77 


12.2 


-2113 




(l)iA 




77 


77 


77 


7.41 


7438 


98 


(3)3n 




77 


77 


77 


17.8 


149 


131 


(4)3S+ 




77 


77 


77 


19.0 


175 


91 






77 


77 


77 


17.9 


157 


134 


(2)in 


D 


Na++Sr(5s4(iiL>) 


-0.52073 


460 


7.32 


3166 


90 


(4)iS+(*) 


E 


77 


77 


77 


8.15 


-2008 


129 


2"^^ well 




57 


77 


77 


14.3 


-1191 


106 


^st barrier 




77 


77 


77 


9.79 


-516 




2"'^barrier 




77 


77 


77 


13.6 


-1158 





TABLE III. Same as Table Ull for NaSr"*". For the (3)^S^ state, the barrier actually refers to the 
position of its pseudocrossing point with the (4)^S^ curve. Due to the peculiar shape of the latter 
potential curve, two barriers are identified. 
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State 


Label 


Diss, limit 


-E'diss(a.u.) 


6{cm ^)) 


Re 








X 


K++Sr(5s2) 


-0.61548 


-184 


7.86 


4977 


81 


(2)iS+ 


A 


K(4s)+Sr+(5s) 


-0.56487 





13.8 


834 


26 




a 




77 


77 


8.75 


6427 


65 


(i)3n 


b 


K++Sr(5s5p3p) 


-0.54890 


-274 


7.27 


4884 


88 


(2)3s+(*) 


c 




77 


77 


8.25 


+327 


73 


2""^ well 




7: 


77 


77 


16.6 


-1115 


22 


barrier 






77 


77 


11.0 


-1801 




(2)3n 




K++Sr(4s3(i3Z)) 


-0.52991 


343 


8.64 


3042 


60 


{If A 




77 


77 


77 


8.26 


3838 


70 


(3)3S+ 




)7 


77 


77 


11.0 


1959 




(i)in 


B 


K++Sr(5s4(iiD) 


-0.52073 


460 


8.34 


4225 


63 


(3)iS+(*) 


C 


77 


77 


77 


8.53 


-5757 


65 


2^^"^ well 




77 


77 


77 


15.2 


-1714 


23 


barrier 




77 


77 


77 


12.0 


-1603 




(l)iA 




77 


77 


77 


8.18 


5361 


72 


(2)in 


D 


K++Sr(5s5pip) 


-0.51779 


-442 


8.32 


2594 


54 


(4)iS+(*) 


E 


77 






13.1 


1205 


48 



TABLE IV. Same as Table [U for KSr"^. The uje value is not reported for the (3)^S"'" curve, as it 
is not significant due to the peculiar shape of the well induced by the interaction with the (2)'^S"'" 
state. 
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State 


Label 


Diss, limit 


Sdiss(a.u.) 


(5(cm-i)) 


Re 




We 




X 


Rb++Sr(5s2) 


-0.61548 


-184 


8.49 


3638 


55 


(2)iS+ 


A 


Rb(5s)+Sr+(5s) 


-0.55886 





14.0 


848 


21 




a 






77 


9.41 


6121 


48 


(i)'n 


b 


Rb++Sr(5s5p3p) 


-0.54890 


-274 


7.90 


3390 


59 


(2)3s+(*) 


c 






77 


8.83 


+1492 


46 


2"'=' well 




5? 


77 


77 


18.0 


-864 


15 


barrier 






77 


77 


11.2 


+2537 


- 


(2)3n 




Rb++Sr(4s3d3L>) 


-0.52991 


343 


9.32 


2122 


42 


(1)3A 










8.90 


2632 


47 


(3)3S+ 






7? 


77 


11.3 


1529 




(i)in 


B 


Rb++Sr(5s4diL») 


-0.52073 


460 


9.0 


3056 


43 


(3)iS+(*) 


C 




77 


77 


9.15 


-2650 


44 


2"'=' well 




5? 


77 


77 


16.0 


-1135 


17 


barrier 




?; 


77 


77 


13.6 


-941 




(l)iA 






77 


77 


8.77 


3989 


50 


(2)in 


D 


Rb++Sr(5s5pip) 


-0.51779 


-442 


9.18 


1554 


33 


(4)iS+ 


E 




77 


57 


13.7 


925 





TABLE V. Same as Table [U for RbSr^". The uje value is not reported for the (3)^S^" and (4)^S"'" 
curves, as it is not significant, due to their peculiar shape. 
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State 


Label 


Diss, limit 


£'diss(a.u.) 


(5(cm-i)) 


Re 




We 




X 


Cs++Sr(5s2) 


-0.61548 


-184 


9.05 


2714 


44 


(2)iS+ 


A 


Cs(6s)+Sr+(5s) 


-0.54845 





13.9 


1133 


20 




a 






77 


9.90 


7003 


41 


(i)3n 


b 


Cs++Sr(5s5p3p) 


-0.54890 


-274 


8.30 


2738 


49 


(2)3s+(*) 


c 




77 


17 


9.38 


+2293 


37 


2"'^ well 




5? 


77 


77 


22.5 


-62 


5 


barrier 






77 


77 


12.0 


-1133 




(2)3n 




Cs++Sr(4s3d3D) 


-0.52991 


343 


9.90 


1450 


33 


{If A 






77 


77 


9.46 


1779 


37 


(3)3S+ 






77 


77 


12.4 


794 




(i)in 


B 


Cs++Sr(5s4(iiD) 


-0.52073 


460 


9.41 


2358 


33 


(3)iS+(*) 


C 




77 


77 


9.67 


-1874 


35 


2"'^ well 




7? 


77 


77 


17.0 


-784 


13 


barrier 




■>■} 


77 


77 


14.2 


-607 




(l)iA 




■>•} 


77 


77 


9.29 


2987 


41 


(2)in 


D 


Cs++Sr(5s5pip) 


-0.51779 


-442 


9.67 


1277 


30 


(4)iS+ 


E 


11 


77 


11 


14.6 


597 


58 



TABLE VI. Same as Table HIl for CsSr+. The value is not reported for the (3)^S"'' curve, as it 
is not significant, due to the peculiar shape of the well induced by the interaction with the (2)'^$]"'" 
state. 



30 



